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RESUMEN 

Del 1 al 4 de febrero de 2016, se realizo en la Universidad de las Fuerzas 
Armadas ESPE, el XXII Curso Internacional de Estructuras, orientado al estudio de 
los Elementos Finitos, tema conocido a nivel mundial porque permite resolver 
problemas complejos en forma sencilla, con la ayuda de un ordenador. 

En la primera clase, Amjad Aref presento el desarrollo numerico del 
Elemento Finito Triangular, conocido como CST (Constant Strain Triangles) y 
posteriormente Enrique Morales indico el uso del programa LS-DYNA. 

En la Universidad de Fuerzas Armadas ESPE, desde el 2009, Roberto 
Aguiar viene desarrollando el sistema de computacion CEINCI-LAB para el analisis 
estatico y dinamico de estructuras, por lo que en este artfculo en primer lugar se 
presenta el marco teorico, muy bien detallado del elemento finito CST y la 
incorporacion de nuevos programas al CEINCI-LAB. Se complementa el estudio 
con la solution de un problema estatico y otro dinamico de una viga de section 
constante en voladizo. 

Palabras Claves: Elemento finito CST. Tension Plana. Sistema de computacion 

CEINCI-LAB 
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SUMMARY 

From February 1 st - 4 th , 2016, the XXII International Course Structures was 
carried out at Universidad de las Fuerzas Armas ESPE, aimed at studying the 
Finite Element Analysis which is known worldwide because it allows to solve 
complex problems in a simple way with the aid of a computer. 

In the first class, Amjad Aref presents the numerical development of the 
Triangular Finite Element which is known as CST (Constant Strain Triangle) and 
later Enrique Morales explained the use of LS-DYNA program. 

At Universidad de las Fuerzas Armas ESPE, Roberto Aguiar has been 
developing the system of CEINCI-LAB computer for static and dynamic structural 
analysis since 2009. The first part of this paper presents a detailed theoretical 
framework of CST finite elements, and them its application through new programs 
in CEINCI-LAB. Finally, the paper presents a solution to the static and dynamic 
problem associated to a cantilever beam with a constant section. 

Keywords: finite element CST. Flat stress. System of CEINCI-LAB computer. 


1. INTRODUCCION 


El metodo del elemento finito se ha convertido en una herramienta poderosa 
en la solucion de problemas de ingenierfa. Las aplicaciones van desde un analisis 
esfuerzo y deformacion de estructuras hasta el analisis de flujo de calor, magnetico 
y otros problemas. Mediante este metodo a las regiones complejas se las discretiza 
en formas geometricas simples que son los elementos finitos. (Onate, 2013) 

Este artfculo es una introduccion, para quienes desean estudiar los 
elementos finitos. Por este motivo se realiza un analisis estatico y dinamico de una 
viga en voladizo, para que sea comprensible su teorfa y los programas que se van 
a desarrollar en MATLAB. 

Se destaca que para resolver una viga se puede emplear el metodo de las 
diferencias finitas (Aguiar, 2011). Es una opcion de calculo, tanto en elementos 
finitos como en diferencias la solucion de la ecuacion diferencial que gobierna el 
comportamiento de vigas a flexion se transforma en la solucion de un sistema de 
ecuaciones lineales. (Aguiar y Mroginski, 2010). 

En este artfculo se va a utilizar el elemento finito denominado: Triangulo de 
Deformacion Constante (CST), en donde la region bidimensional va a ser dividida 
en triangulos de lados rectos. Se permite que cada vertice o nodo se desplace en el 
sentido X, Y. Por lo tanto, cada nodo tiene dos grados de libertad (gdl). 

Para la programacion se va utilizar el software Matlab. Las coordenadas 
nodales se almacenan en un arreglo bidimensional que representa el numero total 
de nodos y sus respectivos grados de libertad, CG. Para el ensamble directo se 
necesita un vector que contenga los grados de libertad en los nodos de cada 
triangulo, VC. Con la matriz de rigidez del elemento triangular y el VC se ensambla 
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la matriz de rigidez total, Kt. A partir de la matriz de rigidez se procede a realizar el 
analisis estatico y dinamico. 


2. ELEMENTO FINITO TRIANGULAR (CST) 

La formulacion bidimensional de elementos finitos mediante el triangulo de 
deformacion constante es el mas sencillo, representa una pequena porcion de la 
estructura en 2D. Este elemento tiene los desplazamientos, las componentes de 
traccion y los valores de la fuerza de cuerpo distribuida en funcion de la position 
indicada por (x, y). Los nodos estan numerados en sentido anti horario, como se 
presenta en la figura 1 , en que el nudo 1 es el vertice del angulo de 90 grados, al 
numerar de esta forma se garantiza que la matriz de rigidez del elementos sera 
definida positiva. 



2.1. Coordenadas reales y naturales 


Las coordenadas (x,y), del elemento finito de la figura 1 , son mapeadas en 
las coordenadas uyv respectivamente, coordenadas naturales. Aref (201 6) 


= a x 4- a 2 x + a 3 y 

a) 

v(x,y) = a 4 + a s x + a 6 y 

(2) 


Los valores de ai - 06 ,de las ecuaciones 1 y 2, son parametros desconocidos 
que varfan dependiendo del elemento. Para determinar los valores se encuentran 
los desplazamientos en sentido X, en cada uno de los nudos o nodos de los vertices 
del triangulo de la figura 1 , con lo que se tiene. 


Uj = uOq.yi) = a t + a 2 x t + a 3 y x 
u 3 = u(x 2 ,y 2 ) = % + a 2 x 2 + a 3 y 2 
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u 5 = u(x 3 ,y 3 ) = a t + a 2 x 3 + a 3 y 3 


En forma matricial, se tiene: 


Ui 

u 3 

u 5 


i *1 y\ 

i *2 y 2 

-i * 3 y 3 - - 



[A] 


Luego, para encontrar las constantes a t se tiene que invertir la matriz A 


{ a } 3 xl — [^] 3^3 { w }3xrl 


(3) 



(x 2 y 3 - x 3 y 2 ) (x 3 y t - x ± y 3 ) {x x y 2 - x 2 yj 


(4) 


Donde, A e es el area del elemento finito triangular. 


2.2. Funciones de forma 


Las funciones de forma o de interpolacion definen la solucion de un 
problema, para un elemento de un portico piano con tres grados de libertad en cada 
uno de sus nudos, se tienen 6 funciones de forma (Aguiar, 2014). Ahora para el 
elemento triangular con los grados de libertad indicados en la figura 1 se tendran 
tres funciones de forma, las mismas que se hallan a partir de la siguiente ecuacion 
en que se supone que las deformaciones del elemento son lineales. 


Donde A/7, N2, N3, son las funciones de forma de desplazamiento y se 
indican en la tabla 1 . Se destaca que se ha remplazado los valores de a t de la 
ecuacion 3 y 4. 


u(x,y ) = % + a 2 x + a 3 y 


Sustituir con la matriz [4] -1 


u(x,y ) = N 1 (x,y)u 1 + N 2 (x,y)u 3 + N 3 (x,y)u s 


(5) 
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Tabla 2 Determinacion de las Funciones de forma 



Donde Aji 1 son los elementos de la matriz inversa de A, que se encuentran 
en la ecuacion (4). Al reemplazar estos valores se halla en forma explicita las 
funciones de forma indicadas en la tabla 2. 


Tabla 2 Funciones de Forma para el elemento CST 

Funciones de Forma 

Ni(x, y) = (x 2 y 3 - x 3 y 2 ) + (y 2 - y 3 )x + (x 3 - x 2 )y 


N 2 (x,y ) = (x 3 y t - x^) + (y 3 - y x )x + (x 3 - xjy 


N 3 (x,y ) = {x x y 2 - x 2 y J + (y 1 - y 2 )x + (x 2 - x x )y 


Los valores de la funcion de forma se definen como la superficie plana, 
que se muestra en la figura 3 

En particular, N1+N2+N3 representa un piano con altura de 1 en los nodos 
1 , 2, 3, es un paralelo al triangulo. 


Ni + N2 + N3 = 1 


( 6 ) 


En este caso, se ha obtenido las funciones de forma en coordenadas reales. 
Una vez que se tienen las funciones de forma el problema practicamente esta 
resuelto. 
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N1(x,y) 



N2(x.y) 3 




Figura 3 Funciones de Forma para elemento CST 


2.3. Matriz B 

Las deformaciones unitarias estan dados por: 

[Ex>£y>Yxy] V) 


Sx ~ dx 


du dN t dN 2 dN 3 

- = ^r Ui+ ^r U3 + ^7 U5 


dv dN 1 dN 2 dN 3 

s y~ dj~W U2 + ^y U 4 + ^y U6 


du dv dN 1 
^ xy dy + dx dy 


Ui + 


d_Ni 

dy 


Uo + 


d_Ns 

dy 


Ur + 


6 N ± 

dx 


U2 + 


d_N 2 _ 
dx 


U 4 . + 


3 

dx 


En forma matricial, 

{E (x,y)} 3x 1 = 


e* (x,y) 
( x,y ) 
Y X y(^X, y) 


3x1 


-1 

CD 

5 

0 

dN _ 2 

0 

a N 3 

0 


rU^> 

dx 

dx 

dx 


U 2 

0 

dN ± 

0 

dN 2 

0 

dN 3 


u 3 > 

dy 

dy 

dy 

< 

U 4 > 

dN t 

dN ± 

dN 2 

dN ± 

dN 3 

dN ± 


u s 

-dy 

dx 

dy 

dx 

dy 

dx 


\U 6 J 


{e(x,y)} 3xl = [B(x,y)] 3 x 6 {q} 6xl 


(8) 
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1 

2^ e 


(y 2 - y 3 ) 

o 

(x 3 - X 2 ) 


0 

(x 3 ~ x 2 ) 

(y -2 - y 3 ) 


(y 3 - yi) o 

0 (x x - x 3 ) 
(Xi - x 3 ) (y 3 - y 4 ) 


(yi - y 2 ) 
o 

(x 2 - Xx) 


r u i\ 


o 

(x 2 - Xx) 

(yi - y 2 ) 


u 2 


1% 

)u 4 




r 5 

\u 6 ) 


De donde la matriz de compatibilidad de deformaciones que sirve para 
calcular la matriz de rigidez es la siguiente. 


1 

yzs 

0 

y 3 i 

0 

yi2 

0 

B = 2A e 

0 

*32 

0 

*13 

0 

*21 

- *32 

y 23 

*13 

y 3 i 

*21 

yi2- 


Notese que la matriz B esta completamente definida, se puede obtener 
directamente para un elemento finito cualquiera, esto facilita notablemente el calculo 
de la matriz de rigidez del elemento. 

2.4. Matriz de rigidez del elemento CST 

La matriz de rigidez del elemento en forma general viene definida de la 
siguiente manera: 


K = 


/ 


B t D Bdv e = t 


l 


B t D BdA 


( 11 ) 


Normalmente la integracion se realiza utilizando la cuadratura de Gauss 
(Aguiar 2016, Aref 2016) pero para el caso del elemento finito CST, la matriz de 
rigidez del elemento se halla con la siguiente expresion. 

K= A e t[B] T [D][B] 


Donde: K es la matriz de rigidez del elemento finito CST; A e es el area del 
elemento finito; t es el espesor del elemento; B es la matriz de compatibilidad 
presenta en la ecuacion (10); D es la matriz que relaciona el vector de esfuerzos 
con el vector de deformaciones. 

Para el caso piano esta matriz es de tres por tres y para el caso de tension 
plana, se tiene: 
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D = 


E 

1 — v 2 


'1 v 

V 1 
0 0 


0 ' 
0 

1 - v 
2 . 


(13) 


Donde E es el modulo de elasticidad del material; v es el modulo de Poisson. 


3. PROGRAMAS PARA EL ANALISIS ESTATICO 

Para la utilizacion del programa, la numeracion de los nudos, elementos y 
grados de libertad se indican en la figura 3. 



Figura 3 Grados de libertad, Numeracion de nudos y triangulos. 


Notese que primero se han numerado los nudos restringidos, luego los 
restantes nudos de la malla, de abajo hacia arriba. Los elementos se han numerado 
de izquierda a derecha y de abajo hacia arriba. Finalmente se indica que cada nudo 
tiene dos grados de libertad un desplazamiento horizontal y un desplazamiento 
vertical; como se va a resolver un problema estatico los grados de libertad se 
numeran de corrido. 


La matriz de rigidez de la estructura, se halla por ensamblaje directo y para 
ello se necesita determinar la matriz de rigidez de cada elemento finito y los vectores 
de colocacion, que no son nada mas que los grados de libertad de los nudos, 
escritos en el orden en que fue deducida la matriz de rigidez del elemento. 

En la figura 4 se aprecia con dfgitos mas grandes la numeracion de los 
nudos y con dfgitos mas pequenos la numeracion de los nudos del elemento finito, 
teniendo en cuenta que se numeran en forma anti horaria, empezando en el vertice 
de angulo recto. 



Figura 4 Numeracion de nudos en sentido anti horario para cada elemento. 

Con estas acotaciones a continuacion se presentan los nuevos programas 
que han sido incorporados al sistema de computacion CEINCI-LAB para el analisis 
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con elementos finitos CST. 

■ dibujo_elemFinito_triangulo(x,y, L, h) 

Esta funcion sirve para graficar todas las divisiones triangulares. 

x Division en el eje x al elemento. 

y Division en el eje y al elemento. 

L Longitud total, en metros, 

h Altura del elemento. 

■ [K]=Kmiembro_ElemFinito_triangulo (L,e,h,x,y,u,E) 

K Matriz de rigidez del elemento triangular. 

E Modulo de elasticidad, en T/m2. 

u Poisson, 

e Espesor. 

h Altura del elemento. 

x Division en el eje x al elemento. 

y Division en el eje y al elemento. 

■ [CG] = cg_elemFinito_ triangulo(x,y) 

CG Matriz de coordenadas generalizas. 

x Division en el eje x al elemento. 

y Division en el eje y al elemento. 


■ [ VC] = vc_ elemFinito_ triangulo( x , y , CG ) 

VC Matriz que contiene vector de colocacion. 
x Division en el eje x al elemento. 

y Division en el eje y al elemento. 

■ [ Kt] = krigidez_elemFinito_ triangulo( x , y , VC , K) 

Kt Matriz de rigidez de la estructura. 
x Division en el eje x al elemento. 

y Division en el eje y al elemento. 

VC Matriz que contiene vector de colocacion. 

K Matriz de rigidez del elemento triangular. 

■ desplazamien f o(x , y , L, h , q ) 

Esta funcion sirve para graficar los desplazamientos de los grados de 
libertad verticales. 

x Division en el eje x al elemento. 

y Division en el eje y al elemento. 

L Longitud total, en metros, 

h Altura del elemento. 
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q Matriz que contiene los desplazamientos. 

4. PROGRAMA PARA EL ANALISIS DINAMICO 


Para el analisis dinamico es necesario numerar al primero o al ultimo las 
coordenadas principales. En este caso se numeran al comienzo, como se indica en 
la figura 5. Se esta considerando que en sentido vertical se tiene un grado de 
libertad. Con esta numeracion no se puede aplicar la triangularizacion de Gauss 
para hallar la matriz de rigidez condensada a las coordenadas principales. 



Figura 5 Grados de libertad, Numeracion de nudos y triangulos. 


Ahora se debe encontrar la matriz que contiene a los vectores de colocacion 
de cada uno de los elementos, para ello se mantiene la misma notation del 
problema estatico pero para que sea didactico se repite la figura 4 en la figura 6. 



Figura 6 Numeracion de nudos en sentido anti horario para cada elemento. 

Los programas, en su mayor parte son los mismos que para el analisis 
estatico sin embargo se los presenta nuevamente, para que sean agrupados en dos 
carpetas, la una para analisis estatico y la otra para analisis dinamico. 

■ dibujo_elemFinito_triangulo(x,y, L, h) 

Esta funcion sirve para graficar todas las divisiones triangulares. 

x Division en el eje x al elemento. 

y Division en el eje y al elemento. 

L Longitud total, en metros, 

h Altura del elemento. 

■ [K]=Kmiembro_ElemFinito_triangulo (L,e,h,x,y,u,E) 


K 

E 


Matriz de rigidez del elemento triangular. 
Modulo de elasticidad, en T/m2. 
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u Poisson, 

e Espesor. 

h Altura del elemento. 

x Division en el eje x al elemento. 

y Division en el eje y al elemento. 

■ [CG\=cg_elemFinito_triangulo_sismo{x,y) 

CG Matriz de coordenadas generalizas. 

x Division en el eje x al elemento. 

y Division en el eje y al elemento. 

■ [ VC] = vc_ elemFinito_ triangulo( x , y , CG ) 

VC Matriz que contiene vector de colocacion. 
x Division en el eje x al elemento. 

y Division en el eje y al elemento. 

■ [Kt g d I ] = krigidez_elemFinito_ triangulo_ sismo(x , y , VC , K) 

Kt Matriz de rigidez 

x Division en el eje x al elemento. 

y Division en el eje y al elemento. 

VC Matriz que contiene vector de colocacion. 

K Matriz de rigidez del elemento triangular, 

gdl Numero de grados de libertad. 

■ [T,fi,OM]=orc/en_e/fif(KL,M) 

KL Matriz de rigidez lateral 

M Matriz de masas. 

T Periodo de vibracion. 

Fi Modos de vibracion. 

OM Frecuencia 


5. APLICACIONES 

5.1. Viga en voladizo 

Se propone calcular el desplazamiento que genera una fuerza de 5T, 
aplicada en el extremo como se indica en la figura 7, a una viga en voladizo de 
longitud de 3 metros. Mediante el uso del programa de elemento finito triangular y 
comparar el resultado mediante la teoria de estructuras. 

Seccion: Base=40cm Material: E=240000kg/cm 2 

Altura=30cm 
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h = 0.40m 


A continuation se indica el programa que se debe elaborar para resolver el 
problema indicado en la figura 7. Utilizando el sistema de computation CEINCI- 
LAB. 


IDatos 

L=3; llongitud en metro 

h=0.4; fAltura en metro 

e=0.3; lespesor en metro 

x =40; lx division de los nudos verticales 

y =10; !y division de los nudos horizontales 

£=2400000; It/m2 


u=0.2; lpoi3son 
I 


dibu j o_e 1 emFini t o_t r i angul o { x , y F L , h) ; 

[ K ] =Kmi erobroE lemF ini tot riangulo [L,e,h r x f y r u r E) ; 
[CG]=cg_eI emFini to_triangulo {x r y) ; 

[VC ]=vc_elemFinito_t riangulo (x f y f CG) ; 

[ E<t ] = tr i g i de z_e 1 emFin i t o_t r i an gul o (x f y f VC f K) ; 
gdl=2* (x+1) * (y+2) ; 

P=zeros (gdl r l); 

P (gdl)=-5; 
q=Kt\P; 

desplazamiento (x f y r L f h F q) ; 

I Calculo por teoria de e structures 

1=0 . 3*0 . 4 A 3/12;EI=E*I;L=3 . 0; P=5; 
des=P*L A 3/(3*EI) ; 


Figura 8 Lineas de Programacion. 

En la figura 9 se indica el mallado con el cual se ha resuelto el problema y 
en la figura 10 se tiene los desplazamientos verticales. Finalmente, se compara el 
desplazamiento vertical en el voladizo con el que se obtiene con la teoria de 
estructuras, se observa que existe una muy buena aproximacion. 
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0 0.5 1 1.5 2 2.5 3 


Figura 9 Mallado de triangulos. 



Figura 10 Desplazamiento en la ordenada 
Resultados 

Elementos finitos: q = 0.01 1 2 m 
Teona de estructuras: q = 0.01 17 m 


5.2. Propiedades dinamicas de viga en voladizo 


Se propone calcular los periodos de vibracion de la viga en voladizo, 
presentada en la figura 1 1 que tiene una longitud de 3 metros. Mediante el uso del 
programa de elemento finito. Comparar el resultado utilizando el elemento finito Q4 
y mediante la solucion de una viga modelada como un sistema continuo. La carga 
que gravita de 2 T/m sirve para determinar la matriz de masas 

Seccion: Base=40cm Material: E=240000kg/cm 2 

Altura=30cm 
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Figura 11 Viga en voladizo. 

Las coordenadas principales son los desplazamientos verticales y se halla 
la matriz de Masas, asociada a esos grados de libertad obteniendo la masa total y 
dividiendo para el numero de coordenadas principales. A continuation se presenta 
el programa desarrollado utilizando los programas de CEINCI-LAB. 


%Datos 

L=3; % lor. git nd er. metro 

h =0.4; %Altura en metro 

e=0.3; %espesor er. metro 

x=20Q; %x division de los nudos verticales 

y=5; %y division de los nudos horisor.tales 

E=24GGG0Q; %t/m2 
u=0 . 2 ; %poisson 

% 

dibujo_elemFinito_tri angul o (x, y r L r h) ; 

[ K ] =Kmi emb r o_E 1 emFir.it o_t r i angul o (L r e f h. f x f y r u, E) ; 
[CG]=cg_elemFinito_tri angul o_sismo (x, y) 

[ VC ] =vc_e lemFinito_tri angul o (x r y, CG) ; 

[Kt gdl l=k:rigidez_elemFinito_tri angul o_sismo (x, VC F K) ; 

% K lateral 

na=x+l ; r.gl=gdl ; 
kaa=Kt { 1 : na, 1 : na) ; 
k:ab=Kt ( 1 : na f na+1 : ngl I ; 
k:ba=k:ab ' ; 

kbb=Kt {na+1 : ngl , na+1 : ngl) ; 

KL=kraa-kab fc inv ( Jcbb) *kba; 

k- Ar.alisis dir.aiaico 

WT=2 *L ; m=WT/ ( (x+l> *9.3) ; 

M=zeros (l,x+l) ; 

□ [for i=l : x+1 

M(l r i)=m; 

L end 

M=diag (M) 

[T r f i r OM ] =o rden_e i g (KL r M) 


Figura 12 Lfnea de Programacion. 
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Para tener una buena aproximacion en el calculo de los periodos de 
vibracion es necesario trabajar con muchos elementos finitos, lo que no sucede 
cuando se trabaja con el elemento finito Q4. En la figura 14 se presenta la malla en 
que se ha considerado 100 divisiones en sentido X, y 5 divisiones en sentido Y. 

Se obtuvo la matriz de rigidez de rigidez asociada a todos los grados de 
libertad y se condenso a las coordenadas verticales, como se observa en la figura 
12. Los primero cinco perfodos de vibracion se presentan en la tabla 3. 



Figura 13 Mallado de triangulos. 


Tabla 3 Periodos de Vibracion de los cinco primeros valores. 


Modos 

1 

2 

3 

4 

5 

Periodos de 
Vibracion 
(segundos) 

0.1134 

0.0192 

0.0074 

0.0042 

0.0028 


• Elemento Finito Q4 

En la figura 14 se presenta la discretization del dominio al trabajar con el 
elemento finito Q4 suavizado por flexion. Notese que se tiene 5 divisiones en X, y 3 
divisiones en Y. Se ha numerado solo las coordenadas verticales, que se considera 
igual a todo lo largo de la Ifnea vertical. (Aguiar, 2016) 



Figura 14 Coordenadas Principales 
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Las coordenadas secundarias son los desplazamientos en cada uno de los 
puntos de la malla, de tal manera que se tienen 25 grados de libertad, se halla la 
matriz de rigidez y se condensa a 5 por 5. La matriz de masas es de5 por 5. 

La numeracion de los elementos se muestra en la figura 15 y en la tabla 4 
se indican los perfodos de vibracion. 



Tabla 4 Periodo de Vibracion mediante Q4 


Modos 

1 

2 

3 

4 

5 

Periodo de 
Vibracion 
(segundos) 

0.1385 

0.0219 

0.0078 

0.0042 

0.0028 


El periodo de vibracion encontrado con el elemento finito Q4 es 0.1385 
segundos y el que se encuentra con el elemento triangular CST es 0.1 134 seg. Se 
trabajo con pocas divisiones con el elemento finito Q4 para entender mejor la 
solucion del problema. 

• Sistema Continuo 

La viga en voladizo de flexion se puede modelar mediante un sistema 
continuo descrito con detenimiento en Aguiar (2012). Ahi se aprecia que la solucion 
de la ecuacion diferencial en derivadas parciales conduce a la siguiente ecuacion. 

1 + cos p Coshp = 0 (14) 

p = aL (15) 

, El A 
Wn 2 = — a 4 

m 


(16) 


Triangulo de Deformacion Unitaria Constante (CST) 
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T = 


2 n 


(17) 


Donde L es la longitud de la viga en voladizo; E es el modulo de elasticidad 
del material; / es el momento de inercia de la viga; m es la masa por unidad de 
longitud. Las cinco primeras rafces de la ecuacion (14) son: 


pi=1 .875; p 2 =4.694; p 3 =7.854; p 4 =1 0.996; p 5 =1 4.1 37 (1 8) 


Para calcular el periodo fundamental asociado al primer modo se tiene: 


= 1 . 875 ; a = 0 . 625 ; W n = 53 . 5826 ; T = 0.1173 seg 

En la tabla 5 se indican los periodos para los cinco primeros modos de 
vibracion. 


Tabla 5 Periodo de Vibracion mediante un sistema continuo. 


Modo 

1 

2 

3 

4 

5 

Periodo de 
Vibracion 
(segundos) 

0.1173 

0.0187 

0.0067 

0.0034 

0.0021 


El periodo fundamental de vibracion hallado con el elemento triangular CST 
es 0.1134 seg., que es muy parecido a la solucion exacta hallada con el sistema 
continuo de 0.1 1 73 seg. 


6. CONCLUSIONES 

Se ha presentado el marco teorico que conduce al calculo de la matriz de 
rigidez de un elemento finito triangular CST, siguiendo el desarrollo numerico 
presentado en Aref (2016). Por otra parte se han desarrollado varios programas de 
computacion para el analisis estatico y dinamico de vigas en voladizo con elementos 
triangulares y se resolvieron dos ejemplos. Del estudio realizado, se desprende: 

■ La contribution al sistema de computacion CEINCI-LAB para la modelacion de 
una viga en voladizo mediante el triangulo de deformacion unitaria constante 
(CST) dan resultados confiables para manipular pero se tiene que tener mucho 
cuidado debido que la malla que se va a generar debe ser muy densa para que 
estos resultados sean buenos, este principio se aplica al analisis estatico como 
al dinamico. Al realizar la modelacion con una malla poco densa los resultados 
son distorsionados y por tanto no se debe utilizar para un posterior calculo. 
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